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ABSTRACT 



In one method of solution to the problem of making up a schedule 

of work in a "job-shop" type industry there arises a syrometric matrix 

the conflicts among the jobs to be done. This thesis is primarily 
concerned with the transformation of this "conflict" matrix to another 

matrix a diagonalized symmetric matrix. That is a matrix having its 

nonzero elements bunched or clustered near the principal diagonal , 

The transformation is considered as the product of simple trsmsformations 
each of v^ich interchange rows or columns of the undiagonalized or 
partially diagonalized matrix, 

A numerical technique is devised for the systemmatic reduction 
of an undiagonalized symmetric mtrix to successively better diagonalized 
natrices . The technique is programmed for use on the NCR 102A Digital 
Computer. 
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CHAPTER I 



INTRODUCTION TO THE PROBLEM UNDER CONSIDERATION 

1 , Introduction 

In a "Job-shop" type industry there are many jobs— job orders 

to be done, each requiring the use of certain of the varied facilities, 
special tools, and specially trained technicians available in the 
industry. Ibe question of how to make up the time schedule of work 
to be done arises. That is which jobs should be time scheduled first 
considering the complexity and use of varied facilities, etc., in the 
industry, Ihen having time scheduled the first one, which job should 
be time scheduled next, and so forth. 

The mathematical problem to be considered in this thesis will assist 
in the determination of an orderly procedure for time scheduling the jobs 
to be done and will provide the order in which jobs should be considered 
for time scheduling. 

2. The Mathematical Problem 

Given a symmetric matrix A, consisting of zero and nonzero elements; 
find a symmetric matrix C whose elements satisfy the following; 

i) The elements of a given column of A become the elements of 
some column of C . 

ii) The elements of a given row of A become the elements of some 
rows of C. 

lii) The nonzero elements of the matrix C are located along the 
principal diagonal of the matrix. That is to say there is a clustering 
of elements near the principal diagonal. 
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Iv) In addition be able to enumerate the correspondence of columns 
of the matrix A vrith columns of the mtrix C. 

The conditions i and ii together are equivalent to the 
following I 

i' ) Columns may be rearranged, maintaining the order of elements, 

il') Rows may be rearranged, maintaining the order of elements. 

For consider any element of the matrix say (i, j). The column i is 
shifted to say the h-th column not necessarily retaining the order of 
elements . Also row j is shifted to say the k-th row not necessarily 
retaining the order of elements, the element that was in position 
(i, J) must now be in position (h, k). However, the same result is 
obtained by shifting the 1-th column to the h-th column maintaining the 
order of elements, and then shifting the j-th row to the k-th row main- 
taining the order of elements. This must be true for all (i,j)» Thus 
i and ii imply i' and ii* , That i' and ii* imply 1 and ii is obvious. 

An assumption is made that the number of nonzero elements in any 
row or column is small compared to the order of the matrix, say one tenth 
the elements are nonzero. If this is not the case the idea of "diagonal- 
ization" or clustering the elements along the principal diagonal has little 
significance . 

The transformation or sequence of transformation used to find the 
matrix C should be applicable to matrices of large order say 62 x 62 or 
even as large as 500 x 500, Since we are considering large order matrices 
it seems logical that the transformation or sequences of transformation 
should be susceptible to programming on a computer, say the NCR 102A 
Digital Computer. 
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Figure 1 General Form of Matrix 



A symmetric matrix A of order n, an n x n array of zero and non- 
zero elements . 'Hie exact value of the nonzero elements is not important — 
the reasons for this are explained in Chapter IV, 

All nonzero elements P^j are located by coordinates (C^,R^.) 
indicating its position in the matrix array is in the C^^th column 
and the Rj'th row. All elements not designated and enuaarated are 
assumed zero. 

The columns are designated numbered in the usual 

fashion from left to right. 
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R numbered from the bottom 
D 



The rows are designated 
upward . This unusual numbering Is associated with the numerical technique 
used to "diagonalize" the matrix. 



have the same value, namely, n + 1. 

The distance, D, of an element from the principal diagonal is 



the jth row, 

4, The Measure of Diagonal izat ion 6f a Jfetrlx 

Some criteria must be set up to measure the clustering or bunching 
of the nonzero elements along the principal diagonal to distinguish 
between the matrices . That is to determine which of two matrices is 
more diagonalized. This measure is relative and serves only as an 
index. 

A simple criterion would be to minimize the distance of the element 
farthest from the principal diagonal. 



This will produce a certain amount of diagonalization , However, a 
simple example will show that a more powerful criterion is needed. 
Figures 2a and 2b both have the same Min Ikx D, namely 2, however. 
Figure 2b is more compact along the principal diagonal. 




1 Sjj -(n * 1) 



= D. Tliis is an integer distance along the ith column or 



Min ffex 



(n . 1) -Sjj 



Min Max D 



|1 0 1 0 0 

10 1 0 1 0 

11 0 1 1 0 
!0 1 1 1 1 
!0 0 0 1 1 



!i 1 0 0 0 
11 1 0 1 0 
,00110 
0 1111 
;_o 0 0 1 1 



(a) 

Figure 2 



(b) 

Examples of Min Max D Criteria 
4 
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Consider the criterion that the sum of the distance from each 



element to the principal diagona] squared be minimized, 
mn / Q n 1) - SjJ — mn 

Each element far from the principal diagonal will contribute greatly 
to the summation, whereas elements next to the principal diagonal will 
have negligible effect on the summation. Thus the greater the cluster- 
ing along the principal diagonal, the smaller will be the summation. 

Ihe previously considered Min Max D criterion will tend to be satisfied 
also . For consider Figure 3 showing two columns (all other columns are 
assumed the same in both cases and , therefore , contribute the same to 
the summation in both cases), 

) 

K 

1 

Figure 3 Example of / for Row Interchange 

In the summation note that only the contribution of the outer element 
in each column is affected by the interchange of the columns under con- 
sideration in Figure 3. In (a) the outer elements contribute to the 
2 2 2 

summation (k+1) + (k-l) or 2k '^♦2. In (b) these same elements con- 

tribute 2k . Thus (b) is better diagonalized in the sense of both 
Min ^fex D and Min > D^. 
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Ckl 
1 0 
1 1 
0 0 



1 1 
1 1 
0 0 



j 

K 






1 0 
0 1 
1 1 



1 1 
1 1 
0 0 



(a) 



(b) 
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For the remainder of the discussion and in testing the merits of 



a transformation under consideration we will use Min 
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CHAPTER II 



SOME FEATURES OF MATRICES .\ND ELEMENTARY TRANSFORMATIONS 
A square matrix A whose elements satisfy the conditions 

possesses symmetry about its principal diagonal, and is called a 
symmetric mtrix. 

If C = B A, then the matrix A is said to be linearly trans- 
formed into the matrix C. The matrix B, which affects the transformation 
of A, may sometimes be advantageously thought of as the multiple product 
of very simple matrices. 



B = B B , 
n n-1 



•® 2®1 



The transformation of A is thus accomplished through a succession of 
simpler transformations, the first of which is the second B^, and 
so forth. 

The matrix can be transformed through postmultiplication as well as 
by premultiplication, i .e . , C = A B' , In which case the matrix 
affecting the transformation is decomposed into components 

I I • 

B' = B, B^ B 

“ 2 n 

t ' 

Where B^ indicates the first simple transformation, B^ the second, and 
so forth. 

It is of interest to determine the simplest fundamental forms into 
which can arbitrary nonsingular matrix B may be decomposed, and to 
interpret the individual transformations which they effect upon the form 
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of the matrix A. These are called elementary transformations and the 
matrices which produce them are called elementary transformation matrices. 
There are three fundamental types of elementary transformations and 
correspondingly three fundamental types of transformation matrices. 

The first type amounts to an interchange of two rows or columns of the 
matrix A. The second type is the addition of the elements of a row or 
column to the corresronding elements of another row or column. The third 
type is the multiplication of any row or column of A by an arbitrary 
nonzero factor. For our purpose it is necessaxy to develop only type 
one to a greater extent. 

The transformation may be affected by multiplying A by a transforma- 
tion matrix. If the desired transformation is intended to affeqt the 
ROWS, A is PREmultiplied by the transformation matrix; if the COLUlWS 
are to be affected, A is POSTnlultiplied by the transformation matrix. 

The elementary transformation matrix is formed from the unit matrix 
by performing on it the same elementary transformation which the desired 
elementary transformation matrix is intended to effect in the matrix A by 
means of (pre-)postmultiplication . Thus I (the unit matrix) with its 
pth and qth rows (columns' interchanged yields a transformation matrix 
which by means of pre- (post- 'multiplication interchanges the pth and qth 
rows (columns) of the matrix A. Apparently the transformation matrix 
has two forms depending on whether it is intended to effect a transforma- 
tion of the rows or columns of A, however, it is to be noted that one 
form is the transpose of the other, and further for type one the two forms 
are identical. We will denote this elementary transformation matrix as 

T where the subscripts indicate the rows (columns) being interchanged, 
p-q 
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As an example, for a fourth order matrix to interchange the second and 
third rows (columns\ 



(1 0 0 0 ’ 

_iooio| ^ „ _ « 

2-3 0 1 0 o' ■ 2-3 23 32 

Rows ijO 0 0 ll Columns 



Consider a symmetric matrix A, pre- and postmultiplied by the 
elementary matrix . The resulting matrix 



= '^i-j" ’^i-J 



is a symmetric matrix with both the ith and jth rows and columns inter- 
changed . Continuing this process we may by the elementary matrix 
Tk-jjj produce the matrix 



A’ ’ = W Vm = Vm A ^i-j Vm 



Thus by a multiple product of elementary transformtion matrices pre- 
and postmultiplying the matrix A, the matrix A is linearly transformed 
into a symmetric matrix C. 

Hie matrix C which may be obtained from the matrix A bv means of a 
finite number of elementary transformations is said to be equivalent to A, 
The equivalence of matrices is a mutual relationship since the trans- 
formations are nonsingular and hence reversible. Thus if C can be 
obtained from A by a succession of elementary transformations, it 
follows that A can be obtained by means of elementary transformations 
from C. The state of equivalence of the two matrices A and C may be 
stated 
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P A Q = C 



where P and Q are nonsingular matrices . 



'‘'ow if A is a symmetric matrix ani Q = where each 

is an elementary tremsformation matrix T, if further 

P = Q^= = \ ’A 

then C is also a symimetric matrix. Since 
C = A Q 

The transpose of C is then 

C = Q A (Q ) = Q A Q. 

A matrix which is its own transpose is symmetric. So C is symmetric. 



The preceding discussion has covered the elementary transformation 
matrices in general . Now to apply these concepts in the notation to be 
used in the numerical technique to be developed in Chapter III. 

In order to have a symmetric matrix after premultiplication by 
an elementary transformation matrix which interchanged the ith and 

jth columns we must multiply by an appropriate row interchanging elemen- 
tary matrix transformation. Since the rows are designated from the 
bottom up, we shall designate the appropriate postraultiplying elementary 

matrix as T — -. Where 
i J 

i-^i = (n + ]) 

J ^ J = (n + 1) 
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Examination of Figure 1 will show that this indeed is the appropriate 
row interchanging matrix. 



Thus A) T~. 

As an example consider the matrix A interchanging first the 2nd and 
5th columns , then the corresponding rows , 





1 0 1 0 l1 


1 


1 1 1 0 oi 


1 


1110 0 


A = 


0 110 1 
11110 


A» = 


0 1 1 0 li 

1 0 1 1 l| 


1 


110 0 1 
10 111 




0 0 110 




0 0 1 1 Oj 


0 0 110 




110 0 1 




1 1 0 0 ij 


1 


0 110 1 



then i = 6 - 2 = 4 
j = 6 - 5 = 1 

Aj = A- Ty = A. 

Aj = ( A ) Ty 

Similarly it can be seen that if the postmultiplication elementary 
transformation matrix is pierformed first then followed by the appropriate 
premultiplication elementary transformation matrix we have 



A 



2 




( A T 



ij 



) 



As an example consider the matrix A interchanging first the 2nd and 
5th rows , then the corresponding columns . 



A = 


10 10 1 
0 110 1 
'l 1 1 1 0 


A' = 


0 0 110 
0 110 1 
11110 


= 


10 10 0 
0 110 1 
11110 


1 


0 0 110 
110 0 1 




10 10 1 

110 0 1 

— J 


2 


0 0 111 
0 10 11 
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4' = 4 Tj5 

^2 = ^41 = ^14 

Aj = ( 4 T,, ) 



Note that at^plied as a premultiplication matrix and as a post- 
multiplication matrix produces different transformed matrices. That is, 



^41 ( '’'25 ’ < ^25* > ^41 

Since we are discussing transformations of symmetric matrices only, 
necessarily these elementary matrix transformations will occur in groups 
of two. So for convenience define a pair of elementary matrix transforma- 
tions as the transformations applied to the rows (columns) and the 
corresponding T7J applied to the columns (rows\ 

The given matrix A has a certain D , say ^ ^o’ 

previous example was ^2. If we apply the transformations 

t = <’'ij '> ’'15 

'C- 2 

the transformed matrix will have a new ^ D ', which for the previous 
example is 3^ . If on the other hand we apply the transformations 



V’lJ <' 



_ ^ 

a different transformed natrix is obtained with a certain ^ D'^ which 
for the previous example is 32. Thus it can be seen that different pairs 
of elementary transformation matrices will produce transformed matrices 

tt II 

diagonalized to a greater or lesser degree. Some pairs of elementary 
matrix transformations may even produce a transformed matrix which is 
less diagonalized in the sense of Min D . 
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CHAPTER III 



PROCEDURE 



1 . General 

Basically the procedure is to apply a pair of elementary trans- 
formation matrices to the undiagonalized or partially diagonalized matrix 



a. If the transformed matrix is better diagonalized, then use this 
matrix as the matrix to be diagonalized and repeat the process applying 
elementary transformation matrices to this partially diagonalized matrix. 

b. If the transformed matrix is NOT better diagonalized, apply a 

different pair of elementary transformation matrices to the undiagonalized 
matrix and test for again . 

We can choose the pair of elementary transformation na trices either 
by some random process or systematically. This thesis is concerned with 
a systematic procedure for improving the diagonallzation of a given 
matrix. 

2, A systematic Procedure 

An element, let us designate it Pjj* which lies farthest from the 
principal diagonal individually contributes the most to 21 D^ of any 
element Pj^j in the matrix. Since by symmetry there are two such elements 
let us consider the one with the largest S, that is Sjj, It appears that 
a pair of elementary matrix transformations which involve this element 
will be a good transformation to hypothesize naking. Furthermore, to move 
this element to a position one unit from the principal diagonal will 
reduce the contribution to of* this element to a minimum. However, 

other elements of the rows and columns interchanged will change distance 



and test the diagonallzation in the sense of minimizing 
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from the principal diagonal so that one must test this hypothesized trsms- 
formation to see if it improves the diagonalization. If it does we make 
the hypothesized transformation and repeat the process again choosing 
the new Pjj with a meiximum Sjj, 

If the previously hypothesized transformation does not improve the 
diagonalization we hypothesize another transformation. If the previous 
elementary transformation pair T^j was a cremultiplication first, i.e., 
row interchange first, we now hypothesize T , making a column interchange 
first such that Pjj is moved to one unit from the principal diagonal. If 
this results in a better diagonalized matrix we make this hypothesized 
transformation and repeat the process again, choosing the new Pjj with a 
maximum Sjj. 

If the hypothesized transformation moving Pjj to a position one 
unit from the principal diagonal does not Improve the diagonalization of 
the matrix, hypothesize moving P^j to two units from the principal 
diagonal. If neither of these pix>duce a better diagonalized matrix, 
hypothesize moving P^j to three units from the principal diagonal and 
continue until the interchange is with the row (column) adjacent to the 
element Pjj. 

If all these hypothesized transformations vdth this element Pjj fail 
to produce a better diagonalized matrix disregard this element and choose 
from among the remaining elements the element with the largest S^j , Repeat 
the process of hypothesizing transformations with this element until one 
is found which improves the diagonalization. If none is found among all the 
elements with greater than n + 1 , i.e., we have looked at all possible 
moves of all ele ents on one half of the symmetric matrix, then the matrix 



U 



which we are attempting to diagonalize is itself the best possible 
diagonalization in the sense of minimizing ^ . 

By hypothesizing moving this "wildest" element to one unit from 
the principal diagonal first, we are maximizing the possible Improvement 
in on-r stage of diagonalization. The actual improvement achieved in one 
stage will usually be somewhat less however. The distance between the 
rows (columns) being interchanged must be variable since it can be shown 
that a system which hypothesizes making the interchange always with an 
adjacent row (column) will not achieve maximum diagonalization. 

Figure 4a is the best diagonalization possible using a hypothesis rule 
of always making the interchanges with the adjacent column (row). 

Figure 4b is the best possible diagonalization using the systematic 
procedure outlined above on the same matrix. 



1010000000 
0100001000 
1011110000 
0011110000 
0011110000 
0011110100 
0100001010 
0000010111 
0000001111 
,_0 0 0 0 0 0 0 1 1 1 _ 

a 



1100000000 
1111100000 
0111100000 
0111100000 
0111110000 
0000111100 
0000011100 
0 000011110 
0000000111 
poooooooii 

b 



Figure 4 Example of Systemmatic Procedure For 
Better Diagonalization 



3, Systematic Procedure Applied to the NCR 102A Computer 

All column and row numbering is done in the octal number system 
for the simple reason that the NCR 102A computes in binary number 
system using data inserted in the octal number system. This 
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OVERALL FLOW CHART-DIAGONALIZATION PROCESS 




Hypothesize interchange moving 
columns first, move P^-j so that 
Y - Djj initi ally 



! i Save these P^ . ' 
as a partially 
Diagonalized 
Matrix 



ZJ 



Return part- 



new \ 

• Yes^ — ^No~ ially ordered 



?ld 



P^j’s to memory 



No- 







Yes 



j Alter Comnands so P^j under consider-' 
- ation is in cell N-1 i 







Figure 5 
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information could be inserted as decimal numbers. However, before 
the computer can start to compute a subroutine would have to be 
employed to convert all data to octal. 

The form and amount of information about a particular nonzero 
element of the matrix A which will be placed in the computer Is dictated 
in part ly the computer and in part a desire by the writer to contain in one 
computer storage word all information relative to this particular element. 
Since a computer storage word is composed of 12 octal digits it is conven- 
ient to store the information about a particular as follows; 
i) for mtrices of order 63 x 63 or less 



l£-?_ 


' i ! 


0 i 


< 


'V'’ ^ '' 





Co «. 



JO 






Since numbers as large as 63 can be represented in 2 octal numbers (63 
decimal is 77 octal) this representation is possible. Where 

Column designation of the nonzero position in the matrix before 
any diagonalization has been effected. 

Row designation of the nonzero position P^j in the mtrix before 
any diagonalization has been effected . 

Current column designation of this P^^j is the matrix at any stage 
of diagonalization . 

Rj Current row designation of this P^^ in the matrix at any stage 
of diagonalization. 

Sj^ j The current C^ Rj . 

It is to be noted that the , R ^ , and S^j may change from time to time 
during the diagonalization process. However, the C^^ and Rj^ remain 
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unchanged and thus a one to one correspondence between the rows (columns) 
of the original matrix and the current partially or completely diagonal- 
ized matrix is obtainable by comparing the and Rj and C^). 

The detailed numerical techniques described later are based on 
stored word structures for order 63 x 63 or less , 



ii) for matrices of order 511 x 511 or less 




Numbers as large as 511 can be represented in 3 octal numbers (511 decimal) 
is 777 octal). However, certain information must be sacrificed in order 
to contain the essential information in one computer storage word length. 
For these larger matrices the value of must be formed vdierever needed 
rather than being immediately available. The s 3 nnbols have the S6ime 
meaning as above. 



Determination of the Largest Sjj. 

In order to determine the "wildest" element or more precisely the 
element with the largest , it is necessary to order the elements 
or at least order some of them, say the largest 6, The reason that more 
them one largest must be obtained is that it may happen that none of 
the hypothesized transformations using the Pjj element with the largest 
Sjj will improve the diagonal Izat ion. However, it is unlikely that 
among the elements with the 6 largest Sj^j's at least one will not have a 
transformation which will improve the dlagonallzatlon. If one of them 
does improve the dlagoned. izat ion the resulting will be reordered 



IS 



FLOW CHART - ORDERING S's 




. Store separately the partially 
ordered P's 
j Initialize X to N 
jTojiext segmen^ of progr am 

Figure 6 
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and the 6 largest of the new matrix obtained. It is to be noted 

that the technique of ordering tends to order not only the 6 largest but 
to relatively orier all the elements . Thus after several partially 
diagonalized matrices are obtained, all the will be relatively well 

ordered . 



The technique used is simply to compare the S^j of the K and K 1 
storage cells, putting the with the smaller in the K storage cell. 
Then compare the of the K + 1 and K + 2 storage cells, putting in 
K + 1 the with the smaller Continuing this process until, the 

N (last) storage cell is reached, then put both away. Decrease 

the number of storage cells under consideration by one, that is, for the 
second time through this ordering process only the first N-1 storage cells 
woul 5 have their ^ ' s compared. To determine if the preassigned number 
of largest S^j's have been obtained, test if N-6 is the last storage cell 
to be considered for the next time through the ordering process . If not 
repeat the process starting with the first two storage cells ; if so go 
on to the next segment of the program. 



faking a Hypothesized Transformation 

For the first hypothesized transformation we consider interchanging 
the elements of two rows first. The P^j stored in the N storage cell will 
have the largest S^j. It is easily determined which row this P^j is in, 
namely the Jth row, and the distance from the principal diagonal, namely D. 
Then J-D is the rov; which has the Ith column on the principal diagonal. 

Let K=J-D+1 

Then the interchange of the Jth and Kth rows will place Pjj next to the 
principal diagonal. Now to maintain s;'mmetry two columns must also be 
interchanged . It is easily seen that the columns 
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Flow CHART - Mh-KING a HYPOTHESIZED TRAIL'S PORNWT ION 




B — 




Find Pjj in Xth ptorii?© cell 

Print P used to determine transformation 
Form Y = S - (n + I) 






Form in buffer: J, K'=J-y*l in last 
2 octal placess 

Form in buffer; i^(n+l^-K, and J=(n+1)-J 
shifted left 2 octal places j 

..lpP'3^P_.^^id_put away for P 1 

• >1 



I Place p" in buf f er 

. JJ-J ) J.nd jK-j ) ^ 



I 

J 






/ \ 



/ \ 

J-j )>0 --»Yes <K-j )>0 /-^Yes 



X 

\ 






r 



Extra ct j=K into P 



No 

1 

iFxtrac t j=J int o T | 



orm 



J-i and K-i 



(j-i)>o'>->Tes 





^(K-i )>0/^es 



No 





fFort-ract i=J into P I 



^ona now S, extract into P 
Form D=vC-.'n+l) 

Form add to partial ^ 

j Put away P 

' ■liter lookup and rut av^-ay addresses 



1 




-Yes^\TM lOy 




-Is F’+l > address of P? 
No 




Yes<— vTM l2>->iJo 

V- 



\ 

Irint 5 ] 






Yes^ 




TM 13 / - >No 

— If previous 
this ^ D2 



'Tally an iteration 
Print iteration and^D^ 



i 1 



Recycle 



I 

I 



A ; 



Figure 7 



Reload memory with 
partially ordered 
P's previously saved. 
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— >Ye: 
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>Yes — 
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No 
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To Correspondence of Original 
and new Columns routine 
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J = (n +1) - J 
K = (n + 1) - K 

are the appropriate columns to interchange to maintain the symmetry. 

These interchanges are the equivalent of the elementary matrix trans- 
formation 

JK ^ JK. 

Forming these quantities J, K, J, K where they are readily available 
proceed to examine each performing the following tests: 
i) if J = J then make J = K 

if J = K then make J = J 
othei*wise row designator remains unchanged, 
ii) if i = J then make i = K 

if i = K then make i = J 

otherwise column designator remains unchanged . 

2 

Determine the new of this element and then D. Form D and add the 
contribution of this element to the partial D already obtained. 

When all the have been so treated, test to determine if the 

V.;" 2 — 2 

D for the hypothesized change is less than ^_^D for the matrix 
before transformation. If indeed this new D is less, we have a 
matrix more diagonalized than before and we start the process over again, 
ordering by the S^^j's, etc. 



If, however, the new D is not less than before discard this 
hypothesized transformation and replace the storage cells with the 
partially ordered set of Pj^j's which had previously been stored separately. 
This time hypothesize a transformation by interchanging elements of two 
columns first. Taking again P^j in the N storage cell determine which 
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column ?jj is in, namely, I, and the distance from the principal diagonal. 
Then I - D is the column which has the Jth row on the principal diagonal. 
Let L = I- D + 1 

Then the interchange of the Ith and Lth columns will place Pjj next to 
the principal diagonal. Now to maintain symmetry two rows must also be 
Interchanged . It is easily seen that the rows 

I = (n + 1) - I 
L = (n + 1) - L 



are the appropriate rows to interchange to maintain the symmetry. These 
interchanges are equivalent to the elementary matrix transformation 
pair Tjj^ A T^. 

Forming these quantities I, L, I» L where they are readily available 
proceed to examine each j performing the following tests : 
i) if J = I then make j = L 

if J = L then make j = I 

otherwise row designator remains unchanged, 
ii) if 1 = I then make i = L 

if i = L t^en make i = I 

otherwise column designator remains unchanged . 

2 

Determine the new of thin element and then D. Form D and add the 
contribution of this element to the partial ^ D^ already obtained. 

When all the P^^ have been so treated test to determine if the 
D for the hypothesized change is lees than ^ D for the matrix 
before transformation. If indeed this new ^ is less, we have a matrix 
more diagonalized than before and we start the process over again, ordering 
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by the etc. 

If, however, the new ^ is not less than the ^ for the matrix 



before transformation, discard the hypothesized transformation, replace 



previously been stored away. 

Now make K' = K 1 
L' = L f 1 

If K' •<C J hypothesize making interchanges as before using K' and L' 
in place of K and L, 

If, however, K' = J, it means that we have hypothesized all possible 
row (colu3Dn) interchanges which shift the Jth row (Ith column) amd none 
of them improve the diagonalization of the matrix. Make N’ = N - 1, 
that is repeat the hypothesizing transformations process using the next 
wildest element. If none of these hypothesized transformations will 
improve the diagonalization of the matrix, make N' = N - 2, and so on, 
continuing this process until all elements have been tried. If none of 
the hypothesized transformations improve the diagonalization of the matrix, 
then the best diagonalization possible using the criteria of minimize 





attained . 



24 



CHAPTER IV 



A PHYSICAL PROBLEM 

Consider a job-shop type industry such as a shipyard which may be 
considered typical . The products may be considered the accomplishment 
of Job orders, each of which may require the use of special facilities, 
special machines and equipment, and the special talents of specialized 
technicians . 

Let us consider several typical or probable job orders: 

Job Order #1: Main Battery Range finder; clean, repaint pris- 

matical marks on lenses , reevacuate rangefinder tube , fill evacuated 
tube with nitrogen, remount and realign the rangefinder. 

Job Order #2! Main Battery Fire Control Equipment; correct 
systemic errors in synchro system between radar and optical portions 
of fire control system, also accomplish electronic modification XX. 

Job Order 0 * Binoculars; collimate and adjust 60 pairs of 
binoculars . 

Job Order #4: Radio Equipment; overhaul model lY receiver. 

And so forth , , , , 

The job orders are analyzed with respect to special machines, 
facilities, and specialized technical work required. For example. Job 
Order §1 will require, l) all the facilities of the optical section for 
work space during cleaning, 2) a specialized technician to repaint the 
prismatical marks on the lenses, 3) special equipment to evacuate the range 
finder tube afterwards, 4) special equipment to refill the range finder 
tube with nitrogen, 5) special alignment equipment smd technical assistance 
during remounting. 
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Job order 
Job order #2 
Job order #3 
Job order #4 
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Job order i^2t while primarily the concern of electronics technicians 
and requiring electronics equipment also will require a synchro expert 
and technicians from the optical shop. 

Job order #3 requires the shop facilities of the optical shop. Thus 
each Job order can be analyzed and a chart prepared as indicated in 
Figure F , where there are m Job orders and n columnar headings , 

Figure 8 in effect represents a m x n matrix of zero and nonzero 
elements. Each x in the figure represents a nonzero element. All the 
nonzero elements of a given row represent the columns (men, machines and 
facili-^ies) which are affected by a particular Job order. Post multipli- 
cation of this m x n matrix M with an n x m matrix, say - the matrix 
transpose- will produce an m xm matrix A. The matrix multiplication is 
carried out in the usual way using the conventions 

(nonzero element) x (nonzero element) = (nonzero element) 

(nonzero element) x ( zero element) = ( zero element). 

This matrix M nF is symmetric with the rows and columns representing 
job order numbers. Thus a nonzero element in the matrix position (i»J) 
indicates that job order i and Job order J have some conflict between them. 

This conflict may be that both Job orders require the use of the space 
facilities of the optical shop or the use of the special nitrogen filling 
equipment. This natrix product MM is thus the "Job order Conflicts matrix. 

This symmetrical conflict matrix is the matrix which is diagonalized 
elsewhere in this thesis . The diagonalization provides an ordering of the 
Job orders such that all job orders which conflict with a given Job order 
are listed in the ordering near that given Job order. Thus a Job order 
which conflicts with no other Job orders will appear at one end of the ordering. 
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Note that the ordering is not unique for the ordering may be enumerated 
in reverse order for the same diagonalized matrix. 

Let us consider again the original job order vs men, machines and 
special facilities matrix M. Since there is nothing special about 
the order that the job orders are listed, let us rearrange the rows 
so that the job orders which conflict with each other are near each 
other. That is precisely the ordering of the job orders given by the 
diagonalization of the "Job Order conflict" matrix. When the matrix 
with the rows ordered in this fashion is scrutinized, it will be noted 
that there is a natural rearrangement of the columns which will yield 
a matrix compact along a line — - for square matrix this line is the 
principal diagonal and for the rectangular matrix this line runs from 
comer to comer. The results will be a matrix such as Figure 9, 

By some criteria such as the type commanders decision that job A 
is most important, the first job order is selected for scheduling. Other 
criteria for finding the starting place are to use the job order row with 
the longest length or the middle row. Having selected the starting row 
by some criteria the complete time scheduling of job order A limits 
certain aspects of Job orders B and C (see Figure 9), The time scheduling 
of B and C in turn limit D and E. This continues until a point is reached 
where F and G are the job orders being time scheduled. They can be 
scheduled independently of each other limited only by the job orders 
already scheduled. This is because they use none of the same men, 
machines or facilities . At this point the remainder of the upper half 
can be scheduled without concern for the lower half causing Interference, 
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APPENDIX 



The complete profi;ram is included here for reference purposes. The 
reader is also referred to the Flow Charts for individual segments of 
the program. 

The program occupies cells 1200 to 1747 in the main memory. The 
program may be operated either minimum access or sequential numbering 
drum. As written the program requires also the use of a magnetic tape 
having blocks 500 to 620 available in both mode 1 6ind mode 2, The nonzero 
elements are stored starting in cell 0000 up to 1077, Channel 11 of the 
memory is reserved for temporary storage during the last section of the 
program— the printout of the correspondence of columns in the undiagonal- 
ized and diagonalized matrices. If part of channel 11 is required for 
the storage of nonzero elements then certain cells and commands must 
be altered for the last segment of the program, (Constant in 1630 t 
and command 16^4 portion) . These could be changed to 1200 vdplng 
out a portion of the program no longer required. If there are less than 
477 (octal) nonzero elements the program can be simply modified to 
eliminate the use of the tape unit. The partially ordered P's then are 
stored starting in 0500, The last iteration is then not saved. 

The only constants which must be filled into program are cells ; 

1440 Number of last cell filled with a nonzero element 

1441 Order of matrix plus one. 
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If P's are initially completely filled out start computing in 1200. 
If only the Co and Ro numbers are in the portion of the P's start 

computing in 1550. If SEN 2020 up program will type out all P's after an 
iteration is completed . SEN 2040 must be up for program to recycle after 
completing an iteration. 
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0000 


n+1 


Order of matrix +1 


N+1 


0000 


0000 


Comparlaon constant 


N 


0000 


0000 


Comparison constant 


0000 


0000 


7777 


Extractor for 


0000 


0000 


7700 


Extractor for C 


0000 


0000 


0050 


Shift control right 8 octal 


0000 


0777 


0000 


Extractor for S 


0000 


0000 


7700 


Extractor for C 


0000 


0000 


0077 


Extractor for R 


0000 


0000 


0001 


Print control} Constant 1 


0000 


0000 


OOlU 


Shift control right U octal 


0001 


0000 


0000 


Print control} Constant 1 


0100 


0000 


0000 


Print control 


5000 


5000 


xxxx 


Comparison constant 


0000 


0000 


Y 


Distance between rows (cols) 
being Interchanged 


0000 


0000 


0050 


Shift control left 8 octal 


2006 


2100 


xxxx 


Comparison constant 


X 


0000 


0000 


Variable constant 


0000 


0000 


0006 


Shift control left 2 octal 


0000 


0000 


0100 




0100 


0000 


0000 


Print control 


0000 


0000 


0006 


Shift control right 2 octal 


0000 


0000 


0010 




7700 


0000 


0000 


Extractor for Co 


7777 


0000 


0000 


Extractor for M^ 


0077 


0000 


0000 


Extractor for Ro 


0000 


0000 


0000 


Tally of Iterations 
^ D^ of last Iteration 


0000 


0000 


0000 


0000 


0000 


0000 


^ D^ of last Iteration 


0000 


0000 


0000 


D^ of current Iteration 


0000 


0000 


0000 


Z. D® of current Iteration 


0000 


0000 


0000 


Buffer save space 


0000 


0000 


0000 


M t1 ft 


0000 


0000 


0000 


M »l If 


0000 


0000 


0000 


tf If If 


0000 


0000 


0000 


If II 11 


0000 


0000 


0000 


ft tf tf 


0000 


0000 


0000 


ft If ft 


0000 


0000 


0000 


tf tf If 
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REMARKS 



ADD 


INST 


M" 




m3 


1510 


21 


1232 


i469 


0001 


1511 


21 


IU73 


2100 


0003 


1512 


21 


1252 


1469 


0001 


1513 


17 


2020 


3000 


1929 


1514 


Ih 


0000 


0100 


0900 


1515 


32 


2100 


l444 


1916 


1516 


09 


3000 


3000 


mx 


1517 


19 


0000 


0102 


0000 


1920 


39 


1916 


1467 


1916 


1921 


3 ^ 


1496 


1916 


1916 


1922 


17 


2040 


3000 


1207 


1923 


22 


0000 


0000 


0000 


I92U 


3 ^ 


3000 


2100 


1200 


1929 


39 


l440 


1492 


2001 


1926 


32 


2001 


1444 


1927 


1927 


21 


0000 


2100 


xxxx 


1930 


21 


1232 


1469 


0001 


1931 


3 h 


3000 


2100 


1914 


1932 


00 


3476 


4179 


6360 


1933 


00 


9664 


6044 


6447 


193*» 


00 


6l4i 


7377 


6641 


1539 


00 


7160 


9949 


4731 


1936 


39 


1476 


1492 


1476 


1937 


39 


2100 


1477 


1477 


19*^ 


3*^ 


3000 


2100 


1979 


19*H 


00 


0000 


0000 


0000 


19*^2 


02 


0000 


0000 


0001 


19»^3 


00 


0000 


0000 


0000 


19*^*^ 


00 


1633 


0000 


0000 


19l»9 


00 


0000 


0000 


0000 


191^6 


00 


0000 


0000 


0000 




00 


0000 


0000 


0000 



Carriage return 
Print Iteration and 
Carriage return 
SEN 2020 UP, print P’s 
Tape block search BL500 
Initialise buffer cmd 
Bviffer load 
Tape wlte 

Add 10 to buffer cmd 

Have we written enough tape? 

SEN 20U0 UP, recycle 

Halt 

After halt, will recycle on 
compute command 
Form N+1 

Extract (N+ll Into print cmd 
Print P’s 
Carriage return 

Flexo code 

•t ft 

II It 

II II 

Tally an overflow 
Wipe out overflow bit 



Buffer save space>mlnaccess 
Shift control right 1 binary 
Buffer save space -mlnaccess 
Constant 

Buffer save space-mlnaccess 
Buffer save space-mlnaccess 



Only the portion, i.e., Co and Ro, need be initially put in 
the computer. If only the portion is filled, starting computing 
in cell 1550 will complete the other portions, namely, current R and 
C and S for all elements P. This portion of the program will also 
compute ar initial . If SEN 2010 up thiszC-D will be entered 
and the initially filled out P's saved in tape unit. 
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ADD INST 



REMARKS 



1550 

1551 

1552 
1595 
155 

1-^1555 

1556 

1557 

1560 

1561 

1562 

1563 

I 56 U 

1565 

1566 
1567 

1570 

1571 

1572 

1573 

157 »^ 

1575 

1576 

1577 
I— 1600 

1601 

1602 

1603 

1 604 
— 1605 



1606 



1607 

1610 

1611 

1612 

1613 
161 U 

1615 

1616 

1617 



M® 



35 


2100 


2100 


1477 


30 


l 44 o 


l 460 


1443 


55 


1445 


1454 


1442 


32 


2100 


1471 


1555 


52 


2100 


1444 


1575 


55 


xm 


2100 


2000 


32 


2000 


1470 


2102 


32 


2000 


1472 


2101 


30 


2001 


1446 


2001 


50 


2002 


1446 


2002 


55 


2002 


2001 


2003 


32 


2005 


1444 


2000 


30 


2001 


1465 


2001 


55 


2002 


2001 


2001 


50 


2001 


i 463 


2002 


52 


2002 


1447 


2000 


30 


2001 


1466 


2001 


36 


2001 


l 44 l 


2001 


26 


2001 


2001 


2001 


35 


1477 


2001 


1477 


37 


1477 


5000 


1556 


35 


2000 


2100 


xm 


55 


1555 


1454 


1555 


35 


1575 


1452 


1575 


34 


1442 


1555 


1555 


21 


1476 


2100 


0002 


32 


l 44 o 


1450 


2100 


35 


2000 


1464 


2000 


32 


2000 


1444 


1456 


17 


2010 


5000 


1434 


22 


0000 


0000 


0000 



00 


0000 


0000 


0000 


02 


0000 


0000 


0050 


00 


0000 


0000 


7700 


00 


0001 


0000 


0000 


00 


0000 


0000 


0002 


00 


0000 


0000 


0077 


00 


0000 


0000 


7777 


00 


0000 


0000 


7700 


00 


0000 


0000 


0006 



COMPLETION OF INFORMATION 
FOR P*8; COMPUTE INITIALED® 
Initialize nartlal Z 
Form N In M^ 

Form N+1 In M^ 

Initialize lookup cmd 
Initialize put avay cmd 
Shift P to htiffer 
Extract Co 
Extract Ro 

Shift Ro 8 octal places 
Shift Co 8 octal places 
^ Extract C and R 
1 Into P 

Shift R 2 octal places 
Form C+R 

Shift S 2 octal places 
Extract S Into P 
Shift S right 2 octal 
Form D=S-(n+l) 

Fonn D^ 

Add D® to partial ^ D® 

Test overflow 
Put away P^ 

Alter lookup 

Alter put away 

Have we owrated on all P? 

Print 

Make comparison constant 
for hirffer-tape cmd 

SEN 2010 UP, shift D® and 
store P's. DO NCT HAVE 
SEN 2040 UP 



CONSTANTS OPTIMUM SPACED 

Shift control right 8 octal 

Extractor for C 

Print control} const 1 

Extractor for R 
Extiractor for 
Extractor for C 
Shift control left 2 octal 
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FLOW CH.IRT • CORRESPONDiENCE OF ORIGINAL NEW COLUMNS 




1 



} Ini tialize lookup and p ut away 



Shift P to buffer, extract S 
Form S-(n+l) = D 



Yes^vD>0 




>No 






ter lookup 



■Put P away in temporary storage! 
lAlter nut_ a way^ (a+l ) | 



— Yea^-<^TM20 > ^Ho 



\ / 



Have we tested all P 
i.e., Ntl > lookup address 







iForin e xtractor ^f __C or ig . 



■»! 






Initia]ize Put at^^y to (a+n) 
Form comparison constants 



[initiaiize lookup P^\ put av/ay in 



j — ^ . ^ 

I jShift to buffer, extract 
[shi ft P^ to buffer, extract 




Yes 



Put P^ in 
P^ storage 



.Alter lookup of P^ 
{Alter put away of 



Put P^ in P^ storage 
^Shift P2 to , S2 to^S^ 



Yes<~ 




TM 21 > 




Have we tested all P’s 
(a+n'>lookup 



No 



Put in last storage posTtxon j 
/liter last storage cell to (a nVj.j 




Is last storage cell > (a vl) 









Form extractor of C new 
Initialize lookup 
* Print heading "Corig-Cne\/' 



Shift P to buffer cell 2000 i 

Extract C, Shift left P octal I 

Put in cell 2001 i 

Print 2 words-first 4 characters only' 
Alter lo okup 



Yes'** 




Is (a ■♦'n)^ Address of P? 




I Alter ends marked * 

, Print Headings 
I "Cnew-Corig‘ 

I Print control to 
: print 2007 and 2000 

j Shifted C is put in 

j 2007 ; jSnd o f ffpg^rami 




i Alter ends marked * 
back to original 



5_settin^s 



Figure 10 



3 ? 



ADD INST 



REMARKS 






m2 



M^ 






© 



1620 


00 


0000 


0000 


0006 


Shift control left 2 octal 


1621 


00 


0000 


0777 


0000 


Extractor for S 


1622 


00 


0000 


0000 


0014 


Shift control left 4 octal 


1625 


00 


0000 


0000 


0077 


Extractor for R 


162 U 


02 


0000 


0000 


0006 


Shift control right 2 octal 


1625 


00 


0000 


0000 


0001 




1626 


00 


0000 


0000 


0077 


Extractor for R 


1627 


00 


2525 


2525 


2525 


Test constant 


1650 


00 


1100 


0000 


1100 


Constant for temp storage 


1651 


00 


1736 


0000 


0000 


Constant 


1652 


00 


0000 


0777 


0000 


Extractor for 8 


1655 


00 


3^36 


6252 


6477 


Flexo code 


163 U 


00 


6560 


7327 


3136 


I! If 


1655 


00 


6252 


6466 


4551 


M II 


-^1656 


50 


l44i 


1622 


2007 


Form (n+l) in M® 


1657 


35 


1447 


2100 


2004 


Shift extractor of C to buffer 


16 U 0 


32 


2100 


1471 


1642 


Initialize lookup 


l64l 


32 


1650 


1444 


l646 


Initialize put away 


r-^l642 


35 


xxxx 


2100 


2000 


Shift P to buffer 


16 U 5 


52 


2000 


2004 


2105 


Extract S into 2005 


1644 


56 


2005 


2007 


2001 


Form S-(n+l)«D 


1645 


3** 


2001 


2100 


1650 — 


Is D greater than 0 


1646 


35 


2000 


2100 


XXXX 


If =0, P is a diagonal element 


1647 


35 


1646 


1452 


1646 


put in temp storage 
Alter put away 


1650 


35 


1642 


1454 


i642^ 


Alter lookup 


J — 1651 


3 ^ 


1442 


1642 


1642 


Have we tested all P*s? 


1652 


35 


1470 


2100 


2002 


Shift extractor of Co to 2002 


—^1655 

1654 


35 

56 


1650 

2000 


l44l 

1615 


2000 

2001 


Form 1100+ (n+l) 
Form 1100+ (n-l) 


1655 


32 


2001 


1444 


1675 


Initialize put away P^ 


1656 


35 


2001 


1452 


2001 


Form 1100+n 


1657 


30 


2001 


l460 


2005 


Shift left 8 octal 


1660 


32 


1650 


1444 


l46i 


Make a comparison constant 
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I 









in 




ADD INST 



PEMARKS 






m2 






1661 
— >1662 
1665 
166k 
1665 

1667 

1670 

1671 
r^<^l672 

I I 

I I-1A7U 

1675 

1676 
i 1677 
' 1700 

1701 

1702 
170*i 
170i* 
1705 

1706 

1707 



55 
52 
52 
55 
52 
55 
52 

5U 
55 
55 
55 
5»* 

55 

56 
56 
5^^ 
5»» 
52 

55 
55 
55 
5*^ 
50 

1710 21 
1711 21 



1712 
1715 
’-^1714 
1715 
j 1716 

! 1717 

I 1720 
j 1721 

' — 1722 




1751 

1752 
1755 

1754 

1755 

1756 

1757 
1740 



55 

52 

55 

52 

50 

21 

55 

21 

54 

55 
57 

52 

52 

52 

54 

55 
52 
52 
52 
22 
00 
00 
00 



1650 

2000 

16^0 

1100 

2006 

XXXX 

2007 

2005 

2007 

1666 

1671 

2005 

2006 
1675 

2005 
1675 
5000 
1671 

2006 

2007 

2005 

5000 

2001 

1655 

1252 

1450 

1650 
xxxx 
2000 
2001 
2000 
1714 
1232 
2005 
1627 
1627 

1651 
1220 
1217 
5000 
2100 
1544 
1252 
1255 
0000 
5456 
4551 

6477 



145 U 
1471 
1444 
2100 
2002 
2100 
2002 
2004 
2100 
1454 
1452 
1666 
2100 
1452 

1454 
l46l 
2100 
1444 
2100 
2100 
2100 
2100 
l460 

1452 

1465 

2100 

1471 

2100 

2002 

l460 

1455 
1454 
1465 
1714 

1627 

5000 

1471 

1471 

l444 

2100 

1627 

1471 

1471 

1444 

0000 

6252 

5156 

6560 



2000 
1666 
1671 

2006 

2104 

2007 

2105 

1702 
XXXX ! 

1666 ^ 

1671 
1666 

XXXX 

1675 
2005 ' 
1662 

1707/D__, 

17054- 

XXXX 

2006 

2004 

1672 
2005 «- 
0005 
0001 
2002 
1714 
2000 
2101 
2001 
0002 
1714 
0001 
1714 

1627 
1751 
1710 
1717 

1716 
1655 
1627^ 

1710 

1717 
1716 
0000 
6466 
6252 
7527 



Make a comparison constant 
Initialize lookup 
Initialize put away 
Shift P^ to buffer 
ExtraOt Into 2004 

Shift P^ to buffer 
Extract Into 2005 

Is C®>C^ 

Put away P^ 

Alter lookup P® 

Alter put away P® 

Tested all P's? 

Put away P^ 

Alter put away P^ 

Alter comparlfon constant 
All P's ordered? 

Prepare put away 
>1 



Put P away 



Shift 

Shift 



P2 

S® 



to P^ 
to 



Prepare a comparison constant 

Print headings 

Carriage return 

Shift extractor for C to 2002 

Initialize lookup 

Shift P to buffer 

Extract C 

Shift C left 8 octal 

Print Co - C 

Alter lookup 

Carriage return 

Have all C's been printed? 

Make tally constant 

Is this 2nd time thru routine? 

Change print control 

Change print control 

Change C put away address 

Wipe out overflow bit 
Change print control 
Change print control 
Change C put away euidress 
Final halt 
Plexo code 
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HI 






